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Abstract 

It is known that in classical fluids turbulence typically occurs at high Reynolds numbers. But can tur¬ 
bulence occur at low Reynolds numbers? Here we investigate the transition to turbulence in the classic 
Taylor-Couette system in which the rotating fluids are manufactured ferrofluids with magnetized nanopar¬ 
ticles embedded in liquid carriers. We find that, in the presence of a magnetic field turbulence can occur at 
Reynolds numbers that are at least one order of magnitude smaller than those in conventional fluids. This 
is established by extensive computational ferrohydrodynamics through a detailed bifurcation analysis and 
characterization of behaviors of physical quantities such as the energy, the wave number, and the angular 
momentum through the bifurcations. A striking finding is that, as the magnetic field is increased, the onset 
of turbulence can be determined accurately and reliably. Our results imply that experimental investigation 
of turbulence can be greatly facilitated by using ferrofluids, opening up a new avenue to probe into the 
fundamentals of turbulence and the challenging problem of turbulence control. 
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Turbulence have been subject to intensive studies as a central issue in the modem science. 
In classical Newtonian fluid, it is well known that the occurrence of turbulence depends on the 
Reynolds number 03 El, and typically in high Reynolds number. Here the Reynolds number 
(Re) is a dimensionless quantity defined as Re = lu/v, where l and u are the characteristic 
length scale and a typical velocity of the flow in a particular geometry, respectively, and v is the 
kinematic viscosity. The quantities l and u may differ essentially from the size of the streamline 
body and mean velocity around, respectively. Therefore for different flow systems the critical 
value of Re can be very different and even approaches infinity, for example, in pipe flows 03. 
In the paradigmatic setting of a uniform flow of velocity u flowing past a cylinder of diameter /, 
when Re is of the order of tens, the flow is regular. For Re in the hundreds, von Karman vortex 
street forms behind the cylinder 01, breaking certain symmetries of the system. Fully developed 
turbulence, in which the broken symmetries are restored, occurs at very high Reynolds number, 
typically in the thousands. In situations where turbulence occurs at high Reynolds numbers, its 
study may be challenging, both experimentally where flow systems of enormous size and/or high 
velocity are required and computationally where unconventionally high resolution in the numerical 
integration of the Navier-Stokes equation is needed. It is thus desirable that fluid turbulence can 
emerge in physical flows of relatively low Reynolds numbers. Meanwhile the ‘way’ that turbulence 
arises can be very different! Consider the three-dimensional Navier-Stokes equations, mainly 
four routes were studied so far |5] [6j]. These present special characteristics and states with (i) 
quasi-periodicity (with two frequencies on two tori) and phase locking, (ii) subharmonic (period 
doubling) bifurcations, (iii) three frequencies (on two- and/or three tori), and (iv) intermittent 
noise. 

Thus it is known that fully developed turbulence may occur after a set of bifurcations [jTJ, as 
in a circular Couette flow with rotating inner cylinder. Moreover, transition to fully developed 
turbulence may depend not only on the Reynolds number, but on the particular characteristics of 
the flow evolutions. For example, in Ref. 03, it was observed experimentally that the flow of a 
sufficiently elastic polymer solutions can become irregular even at low velocity, high viscosity 
and in a small tank that corresponds to very small Reynolds number (on the order of unity). This 
flow has all the main features of fully developed turbulence: a broad range of spatial and temporal 
scales. In general the understanding and “turbulence control” flUEl potentially has a large impact 
on society in particular in economical interests. The energy dissipation of turbulent flows is much 
larger than that of laminar ones. Consequently it is more costly to transport fluid through a vessel 
or to propel a vehicle if the flow is turbulent. For example, in oil pipelines the pressures required 
to pump the fluid is typically up to thirty times larger than would be necessary if the flow could 
be held laminar. Thus to keep the flow laminar until higher Reynolds numbers nowadays it is 
common to add polymers into the flow in oil pipelines. Similar one could think to add ferrofluids 
in circular closed systems to keep the flow laminar due to an applied magnetic field. Studying 
transition to turbulence in non-Newton fluid systems continues to be an interesting topic. 

In this paper, we investigate the Taylor-Couette flow lITOil in finite systems (e.g., aspect ratio 
T = 20) where a rotating ferrofluid ifTH is confined by axial end walls, i.e., non-rotating lids, in 
the presence of an external magnetic field. The classic Taylor-Couette flow of non-ferrofluid has 
been a computational lITK [T3 i| and experimental lfl4Ul6l paradigm to investigate a variety of non¬ 
linear and complex dynamical phenomena, including the transition to turbulence at high Reynolds 
numbers IfTTl . The corresponding ferrofluid system we study consists of two independently ro¬ 
tating, concentric cylinders with viscous ferrofluid filled in between, which has embedded within 
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itself artificially dispersed magnetized nanoparticles. In the absence of the external magnetic field, 
the magnetic moments of the nanoparticles are randomly oriented, leading to zero net magnetiza¬ 
tion for the entire fluid. In this case, the magnetized nanoparticles have little effect on the physical 
properties of the fluid such as density and viscosity. However, when a transverse magnetic field is 
applied, the physical properties of the fluid can be significantly modified llTT[fT 8 ]. leading to drastic 
changes in the underlying hydrodynamics. For example, for systems of rotating ferrofluid lITTl . an 
external magnetic field can stabilize regular dynamical states [fl9]422ll and induce dramatic changes 
in the flow topology [19,1201123]. (In fact, the effects of magnetic field are particularly important 
for geophysical flows If24l - l27l0 . In general, the magnetic field can be used effectively as a control 
or bifurcation parameter of the system, whose change can lead to characteristically distinct types 
of hydrodynamical behaviors [U9l422ll . In this regard, transition to turbulence in magnetohydrody- 
namical (MHD) flows with current-driven instabilities of helical fields has been investigated ||28l . 
There is also a large body of literature on MHD dynamics in Taylor-Couette flows |[29[ - l32ll . Ex¬ 
isting works on rotating ferrofluids lfT9]42Tl l33T - [35ll . however, are mostly concerned with steady 
time-independent flows. Time-dependent ferrofluid flows have been investigated only recently but 
in the non-turbulent regime [|36l . 

The interaction between ferrofluid and magnetic field leads to additional terms in the Navier- 
Stokes equation [32123.34]. O ur extensive computations reveals a sequence of bifurcations lead¬ 
ing to time-dependent flow solutions such as standing waves with periodic or quasiperiodic oscil¬ 
lations, and turbulence. Surprisingly, we find that turbulence can occur for Reynolds numbers at 
least one order of magnitude smaller than those required for turbulence to arise in conventional 
fluids. The occurrence of turbulence is ascertained by a bifurcation analysis and by examining 
the characteristics of the physical quantities such as the energy, the wave number, and the angular 
momentum. We also find that the onset of turbulence can be determined accurately, in contrast to 
classical fluid turbulence where such a determination is typically qualitative and involves a high 
degree of uncertainty [21. Our findings have the following implications: 

• Ferrofluids under magnetic field is a new paradigm for investigating turbulence, especially 
experimentally where the study can be greatly facilitated due to the dramatic relaxation in 
the Reynolds-number requirement. 

• The critical magnetic-field strength for the onset of turbulence can be pinned down precisely, 
possibly leading to deeper insights into the physical and dynamical origins of the transition. 

• Turbulence can be controlled externally, e.g., by an external magnetic field. 

Results 

Ferrohydrodynamical equation of motion. Consider a Taylor-Couette system consisting of 
two concentric, independently rotating cylinders with an incompressible, isothermal, homoge¬ 
neous, mono-dispersed ferrofluid of kinematic viscosity v and density p within the annular gap. 
The inner and outer cylinders of radii R\ and R 2 rotate at angular speeds u\ and uj 2 , respectively. 
The top and bottom end-walls lf37l are stationary and are at distance T(R 2 — R i) apart, where T is 
the non-dimensional aspect ratio. The system can be described using a cylindrical polar coordinate 
system (r, 0, z ) with velocity field (u r , ug, u z ) and the corresponding vorticity V x u = (£, 77 , Q. 
In our preliminary study we set the radius ratio of the cylinders and the parameter Y to typical 
values in experiments, e.g., R\/R 2 = 0.5 and Y = 20. An external, homogeneous magnetic field 
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of strength H x is applied in the transverse x-direction (x — r cos 6). The gap-width d — R 2 — R\ 
can be chosen as the length scale and the diffusion time d 2 /u can serve as the time scale. The 
pressure can be normalized by pu 2 /d 2 , and the magnetic field H and the magnetization M by 
p/p^v/d , where /i 0 is the magnetic permeability of free space. We then obtain the following 
non-dimensionalized equations governing the flow dynamics [ 2211381 : 

(d t + u • V)u - V 2 u + Vp = (M ■ V)H + x (M x H), 

V ■ u = 0. (1) 

The boundary conditions on the cylinders are u(ri,6,z) = (0,f?ei,0) and u (r 2 ,d,z) = 
(0, Re 2 ,0), where the inner and outer Reynolds numbers are Re i = co 2 rid/u and Re 2 = 0 J 2 r 2 d/v, 
respectively, and r i = R \ / ( R 2 — R i) and r 2 = R 2 / ( R> — R i) are the non-dimensionalized inner 
and outer cylinder radii, respectively. To be concrete, we fix the Reynolds numbers at Re i = 100 
and Re 2 = —150 so that the rotation ratio j3 = Re 2 /Re i of the cylinders is —3/2. 

Equation (|T[) is to be solved together with an equation that describes the magnetization of the 
ferrofluid. Using the equilibrium magnetization of an unperturbed state where homogeneously 
magnetized ferrofluid is at rest and the mean magnetic moments is orientated in the direction of 
the magnetic field, we obtain M eq = y H. The magnetic susceptibility of the ferrofluid, y, can be 
determined by Langevin’s formula lf39l . The ferrofluids considered correspond to APG933 Il40)l 
with x = 0.9. The near-equilibrium approximations of Niklas [j34i 41J are small ||M — M eq || 
and small relaxation time | V x u|r -C 1, where r is the magnetic relaxation time. Using these 
approximations, Altmeyer et al. If22il obtained the following equation: 

M-M eq = 4(VxuxH/2 + A 2 SH), (2) 

where 

c n = T / (Vx + rp 0 H 2 / 6/i$) (3) 

is the Niklas coefficient, // is the dynamic viscosity, p 0 the vacuum viscosity, <f> is the volume 
fraction of the magnetic material, § is the symmetric component of the velocity gradient tensor, 
and A 2 is the material-dependent transport coefficient lf38ll . We choose A 2 = 2, which corresponds 
to strong particle-particle interaction and chain formation of the ferrofluid l(38ll . 

Using Eq. ([2]), we can eliminate the magnetization from Eq. ([!} to obtain the following ferro- 
hydrodynamical equation of motion lf22l[38 |: 

(d t + u ■ V)u - V 2 u + Vp M (4) 

c 2 

= - y[HV ■ (F + A 2 §H) + H x V x (F + A 2 SH)], 

where F = (Vxu/2)xH and p M is the dynamic pressure incorporating all magnetic terms 
that can be written as gradients. To leading-order approximation, the internal magnetic field in the 
ferrofluid can be regarded as being equal to the externally imposed field ll2D . which is reasonable 
for obtaining dynamical solutions of the magnetically driven fluid motion. Equation (|4]) can then 
be simplified as 

(d t + u V)u - V 2 u + X7 Pm = s 2 x {V 2 u - 2A 2 [V ■ (SH)] (5) 

-H x [V x (V x u x H/2) - H x (V 2 u) +A 2 V x (SH)]} . 
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FIG. 1: Bifurcations with Niklas parameter s x : (a) time-averaged modal kinetic energy, (b) its m = 2 
contribution, and (c) spatiotemporally averaged axial flow field at midgap for 6 = 0 and 0 = tt/2 in respect 
of the applied magnetic field. Open and filled symbols are for steady-state and time-dependent solutions, 
respectively. 

This way, the effect of the magnetic field and the magnetic properties of the ferrofluid on the ve¬ 
locity field can be characterized by a single parameter, the magnetic field or the Niklas parameter, 


2(2 + x)H x cn 

Sx ~ (2 + x) 2 ~ X 2 V 2 ( 

Transition to turbulence in ferrofluids. We first present a sequence of bifurcations with the 
magnetic parameter s x , eventually leading to turbulence. For s x ^ 0, the rotating ferrofluid flow 
is intrinsically three-dimensional [IT914221 with increased complexity as s x is increased from zero. 
Figures [TJa-c) show, respectively, three key quantities versus s x : the time-averaged modal kinetic 
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energy defined as 


__ / p2ir pT/2 /•To \ 

y ^E m = ( I / / u m u^rdrd^d6» ) (7) 

V 0 ^-r/2 irj / t 

its rn — 2 contribution E 2 , and the spatiotemporally averaged axial flow field <“*>L at midgap 
for 6 = 0 (along the magnetic field) and 6 = 7t/ 2 (perpendicular to the field). We see that, for 
s x < 0.784, the flow is time-independent. In particular, for 0 < s x < s^ t ~ 0.642, the flow 
patterns show wavy vortices with a two-fold symmetry 09l . denoted as WVF(2), which appears 
in Fig. [2^a) as a “two-belly” structure. As s x is increased through sjf*, this symmetry is broken 
and a somewhat tilting pattern in the wavy-vortex structure emerges [denoted as WVF(2),J, as 
shown in Fig. |2jb). The intuitive reason is that the magnetic force is downward for fluid flow 
near the annulus and upward where the flow exits, resulting in a split in (u z ) 2 zt for 6 = 0 and 
7 t/2. As s x is increased through the critical point s p x ~ 0.784, the flow becomes time-dependent. 
F° r s x > s p , the flow is time periodic (limit cycle), corresponding to standing waves with axial 
oscillations of characteristic frequency ui\, which is denoted as SWO p . Near the onset of SWO p , 
the kinetic energy values E and E 2 are unaffected, as can be seen from Figs.|T[a) and[7[b), respec¬ 
tively. For s x > s qp ~ 0.81, the periodic solution becomes unstable due to the emergence of the 
second incommensurate frequency u 2 (~ cui/5) associated with defect propagation through the 
annulus in the axial direction, leading to a transition to quasiperiodic flow (denoted as SWO pp ). 
[See movie files moviel.avi, movie2.avi, movie3.avi, and movie4.avi in Supplementary Materials 
(SMs), where the additional frequency can be identified visually in the defect starting from the top 
of the oscillating region, propagating downward toward the bottom, and vanishing there.] While 
E continues to decrease through this transition, both E 2 and (u z ) 2 t reach their respective max¬ 
ima at the transition point. As s x is increased further, the flow becomes more complex. Onset of 
turbulence occurs for s x ~ sfl ~ 0.825, where globally the flow starts to rotate in the azimuthal 
direction, as shown in Fig. [2jd) (see also movies movie5.avi, movie6.avi, and movie7.avi in SMs). 
The observed route to turbulence coincides completely with that established previously 01142] for 
conventional fluid at high Reynolds numbers. 

Wave structure of the flows. We next analyze the wave structure of the flows to probe into 
the transition to turbulence in the low Reynolds-number regime. To get numerical solutions of 
the ferrohydrodynamical system, we solve the ferrohydrodynamical system by combining a finite- 
difference, time-explicit method of second order in (r, z ) with Fourier spectral decomposition in 6 
(See Method). The variables can be written as 

777 max 

f(r,9, z,t) = Y fm(r,z,t)e lTne , (8) 

m= TTlmax 

where / denotes any one of (w r , ug, u z , p}. We then perform an axial Fourier analysis 091 of the 
mode amplitudes of the radial velocity field ( u r ) m (z,t ) at midgap to obtain the behavior of the 
wave number k, as shown in Fig. [3} Prior to the onset of turbulence, the number of vortices in 
the bulk fluid is invariant, regardless of whether the flow is time-independent or time-dependent. 
However, due to the symmetry breaking induced by the magnetic field, the spatial structures of 
the flows can differ significantly. The top panel of Fig. [3] shows that the axial wave number k 
begins to split into two branches at s x = sfl', where k is enhanced and reduced as the field enters 
into and exits perpendicularly from the annulus, respectively, and there are dramatic changes in 
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(a) WVF(2) (b) WVF(2) 

s = 0.5 s - 0.7 

X X 


(c)SWO' (d) turbulence 

s = 0.8 s - 0.9 

X X 



FIG. 2: (a-d) For four values of s x cones ponding to WVF 2 , WVF/, SWO p , and turbulence regimes, respec¬ 
tively, isosurfaces of azimuthal vorticity 77 = d z u r — d r u z and contours of the radial velocity u. r (0. z ) on an 
unrolled cylindrical surface in the annulus at midgap. Red (dark gray) and yellow (light gray) colors denote 
77 = ±100 for isosurfaces, and inflow and outflow for contour plots, respectively. While the pattern in (a,b) 
are stationary, the ones in (c,d) are snapshots due to time dependence of the corresponding flow. 

the flow profiles in which the dynamics tend to focus on the 6 « 0 region. As s x is increased 
further, k in the bulk decreases with nearly constant difference between the cases of 6 = 0 and 
6 = 7r/2. However, the kinetic energies and the axial wave numbers are nearly unaffected when 
the flow becomes time-dependent as s x is increased through s p x . As the quasiperiodic regime is 
reached, the energies and the axial wave numbers begin to change where, as shown in Fig. [3j 
the wave numbers for 6 = 0 and 6 = 7t/2 are minimized at s x about 0.8. For both the periodic 
and quasiperiodic regimes, there is little variation in k and only the vortex pair surrounding the 
oscillatory region expands and shrinks periodically. Especially, in the periodic regime the flow 
profile at 6 = 7t/2 exhibits little dependence on s x except for an increasing amplitude, but in 
the quasiperiodic regime the flow shows a strong dependence on s x , due to the emergence of the 
second frequency. As is passed, the flow starts to rotate and the split in the axial wave numbers 
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(a)s = 0.5 (b) 0.7 


(c) 0.8 (d) 0.815 (e) 0.9 
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FIG. 3: Top panel: variation with s x of the axial wave number k in the directions along (6 = 0) and 
perpendicular to {6 = n /2) the magnetic field, (a-e) Snapshots of axial velocity u z for 6 = 0 (dashed 
lines) and 6 = 7t/2 (solid lines) in the annulus at the midgap location for five different values of s x , where 
(c-e) correspond to periodic, quasiperiodic, and turbulent flows, respectively. The presented patterns are 
for m = 0 so that we can identify the largest contribution in the axial Fourier spectrum of (u z )o(z,t). 
In principle one can also identify k from the axial profiles in the figure that gives the axial wavelength A 
and consequently the wavenumber k = 2tt/X. See also movie tiles movie8.avi, movie9.avi, movie2.avi, 
movie4.avi and movie7.avi in SMs. 
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in the directions parallel and perpendicular to the magnetic field vanishes. In fact, as turbulence 
sets in no distinct wave numbers can be identified, as shown in Figs. [2] and ^e) (see also movie 
file movie7.avi in SMs). 
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FIG. 4: For s x = 0.8 (periodic regime), isosurfaces of the relative angular momentum rug — J ’J" rugdt at 
eight time instants in one period of oscillation, where r « 0.058 and the isolevels are rug = ±5 (see also 
movie files moviel0.avi, moviel l.avi, and moviel2.avi in SMs). 

Behavior of the angular momentum. We now examine the behavior of the angular momentum 
in relation with the flow bifurcation sequence and transition to turbulence. Figure [4] shows, for 
flows in the periodic regime, the isosurfaces of differences in the angular momentum rug between 
the full flow and its long-time averaged value, namely, rug—fj rugdt. Onset of the periodic regime 
can be identified by the occurrence of periodic oscillations (up and down) of a single outward 
directed jet of angular momentum rug (see movie files movie2.avi, and moviel.avi in SMs). From 
Fig. [4] we see that, while the full flow pattern in the periodic regime exhibits little variation (see 
movie files moviel0.avi and moviel l.avi in SMs), there is relatively large variation in the behavior 
of the angular momentum (see moviel2.avi in SMs). For example, there is a downward flow from 
6 = 0 and an upperward flow from 9 — n, with opposite angular-momentum values. In the middle 
of the finite system where the effects of the Ekman boundary layers are minimal, the difference 
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in the angular momentum reaches maximum. As the system enters into the quasiperiodic regime, 
the appearance of the incommensurate frequency cu 2 signifies a kind of defects in the propagation 
pattern of the flow. As the flow begins to rotate in the azimuthal direction, turbulence sets in. 



(c) 


r 


(d) 


r 



FIG. 5: (a) Unnormalized and normalized (by the total kinetic energy) maximum cross-flow energy E r ( !j ax 
versus the Niki as parameter s x , where E c f [r. t ) = (v% + is averaged over the surface of the con¬ 

centric cylinder, (b-d) Spatiotemporal evolutions of E c f(r,t ) for s x = 0.8, 0.82, and 0.9, corresponding to 
periodic, quasiperiodic, and turbulent regimes, respectively. Red (yellow) color indicates high (low) energy 
values. 

Low Reynolds number turbulence. An extremely challenging issue in the study of turbulence 
is precise determination of its onset as a system parameter, e.g., the Reynolds number, is changed. 
For the type of low-Reynolds number turbulence in ferrofluid uncovered in this Letter, this 


10 














can actually be achieved. In particular, the transition to turbulence coincides with the onset of 
azimuthal rotation from the quasiperiodic flow. To illustrate this, we consider the cross-flow 
energy E c f(r,t), a commonly used indicator for the onset of turbulence in the Taylor-Couette 
system ll43l . which is the instantaneous energy associated with the transverse velocity component 
at radial distance r. Small (large) values of E c f indicate laminar (turbulent) flows. We find that, 
the maximum value of E c f(r,t), denoted by Efj ax , assumes near zero values for s x < s r j]’ but 
it increases dramatically as s x is increased through s®\ as shown in Fig. [5} a). In this figure we 
also plot the scaled maximum cross-flow energy value, Efj ax /E, versus s x , and we observe that 
E'"j ax and Ef ax / E exhibit essentially the same behavior as the system passes through the onset 
of turbulence. Examples of the spatiotemporal evolution of E c f(r,t ) are shown in Figs. |5}b,c,d) 
for periodic, quasiperiodic, and turbulent regimes, respectively. We observe regular patterns in 
the former two regimes, but no apparent patterns in the turbulent regime. 

Discussion 

To summarize, we have discovered that in the Taylor-Couette ferrofluid system driven by a 
magnetic field, where flow can exhibit axial oscillations but not rotations in the azimuthal 
direction in the regular regime, turbulence can generically arise and its onset can occur for low 
values of the Reynolds number. This is substantiated by extensive computations of the underlying 
ferrohydrodynamical equation through a systematic bifurcation analysis and characterization 
of behaviors of physical quantities. The implications, besides the surprising phenomenon of 
turbulence at very low Reynolds numbers and consequently facilitation of experimental study of 
turbulence, lie in the perspective of controlled generation of turbulence through variations of the 
external magnetic field, making it possible to locate the onset of turbulence with high precision. 
We expect these findings to have values for experimental as well as theoretical investigation of 
turbulence. 

Method 

Numerical scheme for ferrodynamical equation. System @ can be solved [ 19, [2311 221 by 
combining a second-order finite-difference scheme in (r, z) with Fourier spectral decomposition 
in 6 and (explicit) time splitting. The variables can be written as 

f(r,0,z,t)= fm(r, z, t) e imd , (9) 

m= m max 

where / denotes one of {u r , uq, u z , p}. For the parameter regimes considered in our preliminary 
study, the choice m max = 10 provides adequate accuracy. We use uniform grids with spacing 
5r = 5z = 0.05 and time-steps 5t < 1/3800. For diagnostic purposes, we also evaluate the 
complex mode amplitudes f m ,n{j' 1 1) obtained from the Fourier decomposition in the axial direction 

/ m (r, z, t ) = fm,n{r, t)e mkz , where k is the axial wavenumber. 
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Figure legends 

Figure 1: Bifurcations with Niklas parameter s x : (a) time-averaged modal kinetic energy, (b) 
its rn — 2 contribution, and (c) spatiotemporally averaged axial flow field at midgap for 6 = 0 and 
9 = 7t/ 2 in respect of the applied magnetic field. Open and filled symbols are for steady-state and 
time-dependent solutions, respectively. 

Figure 2: (a-d) For four values of s x corresponding to WVF 2 , WVF,. SWO p , and turbulence 
regimes, respectively, isosurfaces of azimuthal vorticity q = d z u r — d r u z and contours of the radial 
velocity u r (9, z ) on an unrolled cylindrical surface in the annulus at midgap. Red (dark gray) and 
yellow (light gray) colors denote q = ±100 for isosurfaces, and inflow and outflow for contour 
plots, respectively. While the pattern in (a,b) are stationary, the ones in (c,d) are snapshots due to 
time dependence of the corresponding flow. 
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Figure 3: Top panel: variation with s x of the axial wave number k in the directions along 
(9 = 0) and perpendicular to ( 9 = tt/2) the magnetic field, (a-e) Snapshots of axial velocity 
u z for 9 = 0 (dashed lines) and 9 = n/2 (solid lines) in the annulus at the midgap location 
for five different values of s x , where (c-e) correspond to periodic, quasiperiodic, and turbulent 
flows, respectively. The presented patterns are for m — 0 so that we can identify the largest 
contribution in the axial Fourier spectrum of (u z ) 0 (z, t ). In principle one can also identify k from 
the axial profiles in the figure that gives the axial wavelength A and consequently the wavenumber 
k = 2tt/X. See also movie files movie8.avi, movie9.avi, movie2.avi, movie4.avi and movie7.avi 
in SMs. 

Figure 4: For s x = 0.8 (periodic regime), isosurfaces of the relative angular momentum rug — 
fo rugdt at eight time instants in one period of oscillation, where r ~ 0.058 and the isolevels are 
rug = ±5 (see also movie files moviel0.avi, moviell.avi, and moviel2.avi in SMs). 

Figure 5: (a) I In normalized and normalized (by the total kinetic energy) maximum cross-flow 
energy E™j? x versus the Niklas parameter s x , where E c f(r,t ) = (u 2 + u 2 z )a( r ) is averaged over 
the surface of the concentric cylinder, (b-d) Spatiotemporal evolutions of E c f(r,t ) for s x = 0.8, 
0.82, and 0.9, corresponding to periodic, quasiperiodic, and turbulent regimes, respectively. Red 
(yellow) color indicates high (low) energy values. 
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